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Abstract 

In an attempt to describe cork-pulling, we model a cork as an 
incompressible rubber-like material and consider that it is subject to 
a helical shear deformation superimposed onto a shrink fit and a simple 
torsion. It turns out that this deformation field provides an insight 
into the possible appearance of secondary deformation fields for special 
classes of materials. We also find that these latent deformation fields 
are woken up by normal stress differences. We present some explicit 
examples based on the neo-Hookean, the generalized neo-Hookean, 
and the Mooney-Rivlin forms of the strain-energy density. Using the 
simple exact solution found in the neo-Hookean case, we conjecture 
that it is advantageous to accompany the usual vertical axial force 
by a twisting moment, in order to extrude a cork from the neck of a 
bottle efficiently. Then we analyse departures from the neo-Hookean 
behaviour by exact and by asymptotic analyses. In that process we are 
able to give an elegant and analytic example of secondary (or latent) 
deformations in the framework of nonlinear elasticity. 

1 Introduction 

Rubbers and elastomers are highly deformable solids which have the remark- 
able property of preserving their volume through any deformation. This 
permanent isochoricity can be incorporated into the equations of continuum 
mechanics through the concept of an internal constraint, here the constraint 
of incompressibility. Mathematically, the formulation of the constraint of in- 
compressibility has led to the discovery of several exact solutions in isotropic 
finite elasticity, most notably to the controllable or universal solutions of 
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Rivlin and co-workers (see for example Rivlin (1948)). Subsequently, Er- 
icksen (1954) examined the problem of finding all such solutions. He found 
that there are no controllable finite deformations in isotropic compressible 
elasticity, except for homogeneous deformations (Ericksen 1955). The im- 
pact of that result on the theory of nonlinear elasticity was quite important 
and long-lasting, and for many years a palpable pessimism reigned about the 
possibility of finding exact solutions at all for compressible elastic materials. 
Then Currie & Hayes (1981) showed that one could obtain interesting classes 
of exact solutions, beyond the homogeneous universal deformations, if one 
restricted their attention to certain special classes of compressible materials. 
A string of results about the search for exact solutions in nonlinear elastic- 
ity followed. Now a long list exists of classes of exact solutions which are 
universal only relative to some special strain-energy functions (for a recent 
presentation of such classes see Fu and Ogden (2001).) These solutions can 
help us to understand the structure of the theory of nonlinear elasticity and 
to complement the celebrated solutions of Rivlin. 

In the same vein, some recent efforts focused on determining the maximal 
strain energy for which a certain deformation field, fixed a priori, is admis- 
sible. This is a sort of inverse problem: find the elastic materials (that is, 
the functional form of the strain-energy function) for which a given defor- 
mation field is controllable (that is, for which the deformation is a solution 
to the equilibrium equations in the absence of body forces). A classical ex- 
ample illustrating such an approach is obtained by considering deformations 
of anti-plane shear type. Knowles (1977) shows that a non-trivial (non- 
homogeneous) equilibrium state of anti-plane shear is not always (univer- 
sally) admissible, not only for compressible solids (as expected from Erick- 
sen's result) but also for incompressible solids (Horgan (1995) gives a survey 
of anti-plane shear deformations in nonlinear elasticity). Only for a spe- 
cial class of incompressible materials (inclusive of the so-called 'generalized 
neo-Hookean materials') is an anti-plane shear deformation controllable. 

Let us consider for example the case of an elastic material filling the annu- 
lar region between two coaxial cylinders, with the following boundary value 
problem: hold fixed the outer cylinder and pull the inner cylinder by apply- 
ing a tension in the axial direction. It is well established that a solution to 
this problem, valid for every isotropic incompressible elastic solid, is obtained 
by assuming a priori that the deformation field is a pure axial shear. Now 
consider the corresponding problem for non-coaxial cylinders, thereby losing 
the axial symmetry. Then it is clear that we cannot expect the material to 
deform as prescribed by a pure axial shear deformation. Knowles's result 
tells us that now the boundary value problem can be solved with a general 
anti-plane deformation (not axially symmetric) only for a subclass of incom- 
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pressible isotropic elastic materials. Of course this restriction does not mean 
that for a generic material it is not possible to deform the annular material 
as prescribed by our boundary conditions, but rather that in general, these 
lead to a deformation field which is more complex than an anti-plane shear. 
And so, we expect secondary in-plane deformations. 

The theory of non-Newtonian fluid dynamics has generated a substantial 
literature about secondary flows, see for example Fosdick & Serrin (1973). 
In solid mechanics it seems that only Fosdick & Kao (1978) and Mollica 
& Rajagopal (1997) produced some significant and beautiful examples of 
secondary deformations fields for the non-coaxial cylinders problem, although 
this topic is clearly of fundamental importance, not only from a theoretical 
point of view but also for technical applications. 

In this paper we consider a complex deformation field in isotropic incom- 
pressible elasticity, to point out by an explicit example the situations just 
evoked and to elaborate on their possible impact on solid mechanics. Our 
deformation field takes advantage of the radial symmetry and therefore we 
find it convenient to visualize it by considering an elastic cylinder. 

Let us imagine that a corkscrew has been driven through a cork (the 
cylinder) in a bottle. The inside of the bottleneck is the outer rigid cylinder 
and the idealization of the gallery carved out by the corkscrew constitutes the 
inner coaxial rigid cylinder. Our first deformation is purely radial, originated 
from the introduction of the cork into the bottleneck and then completed 
when the corkscrew penetrates the cork (a so-called shrink fit problem, which 
is a source of elastic residual stresses here). We call A, B the respective inner 
and outer radii of the cork in the reference configuration and r\ > A, r2 < B 
their current counterparts. Then we follow with a simple torsion combined 
to a helical shear, in order to model pulling the cork out of the bottleneck in 
the presence of a contact force. Figure [T] sketches this deformation. 

Of course we are aware of the shortcomings of our modelling with re- 
spect to the description of a 'real' cork-pulling problem, because no cork is 
an infinitely long cylinder, nor is a corkscrew perfectly straight. In addition, 
traditional corks made from bark are anisotropic (honeycomb mesoscopic 
structure) and possess the remarkable (and little-known) property of having 
an infinitesimal Poisson ratio equal to zero, see the review article by Gibson 
et al. (1981). However we note that polymer corks have appeared on the 
world wine market; they are made of elastomers, for which incompressible, 
isotropic elasticity seems like a reasonable framework (indeed the documen- 
tation of these synthetic wine stoppers indicates that they lengthen during 
the sealing process). We hope that this study provides a first step toward 
a nonlinear alternative to the linear elasticity testing protocols presented in 
the international standard ISO 9727. We also note that low-cost shock ab- 



3 



Figure 1: Shrink fit of an elastic tube, followed by the combination of simple 
torsion and helical shear. (The figure does not respect scales among the 
various deformations) 

sorbers often consist of a moving metal cylinder, glued to the inner face of 
an elastomeric tube, whose outer face is glued to a fixed metal cylinder (Hill 
1975). 

The plan of the paper is the following. Section 2 is devoted to the deriva- 
tion of the governing equations and to a detailed description of the boundary 
value problem. In §3 we specialize the analysis to the neo-Hookean strain- 
energy density, and find the corresponding exact solution. We use it to show 
that it is advantageous to add a twisting moment to an axial force when ex- 
truding a cork from a bottle. The neo-Hookean strain-energy density is linear 
with respect to the first principal invariant of the Cauchy-Green strain ten- 
sor. It is much used in Finite Elasticity theory, although it captures poorly 
the basic features of rubber behaviour (Saccomandi 2004). We thus inves- 
tigate the consequence of departing from that strain-energy density. First, 
in §4 we consider the generalized neo-Hookean strain-energy density — non- 
linear with respect to the first principal invariant of the Cauchy-Green strain 
tensor — to show that in this case torsion is explicitly present in the solution 
for the axial shear displacement, but it is a second-order dependence. Next 
we consider in §5 the Mooney-Rivlin strain- energy density — linear with re- 
spect to the first and second principal invariants of the Cauchy-Green strain 
tensor — and find that it is then also possible to obtain an exact solution 



to our boundary value problem. Its expression is too cumbersome to ma- 
nipulate and we resort to a small parameter asymptotic expansion from the 
neo-Hookean case. Section 6 concludes the paper with some remarks on the 
limitations of the semi-inverse method. 



2 Basic equations 



Consider a long hollow cylindrical tube, composed of an isotropic incom- 
pressible nonlinearly elastic material. At rest the tube is in the region 

A < R < B, < < 2tt, -oo < Z < oo, (2.1) 

where (R, 0, Z) are the cylindrical coordinates associated with the unde- 
formed configuration, and A and B are the inner and outer radii of the tube, 
respectively. 



2.1 Equilibrium equations 



Consider the general deformation obtained by the combination of radial di- 
latation, helical shear, and torsion, as 



r = r(R), 9 = Q + g(R) + rZ, 



XZ + w(R), (2.2) 



where (r, 8, z) are the cylindrical coordinates in the deformed configuration, 
r is the amount of torsion and A is the stretch ratio in the Z direction. Here, 
g and w are yet unknown functions of R only (The classical case of pure 
torsion corresponds to w = g = 0, see the textbooks by Ogden (1997) or by 
Atkin and Fox (2007), for instance.) 

Hidden inside (12.21) is the shrink fit deformation 



r(R), 



6 = Q, 



XZ, 



(2.3) 



which is (I2.2p without any torsion nor helical shear (r = g = w = 0). 

The physical components of the deformation gradient F and of its inverse 
F^ 1 are then 
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respectively. Note that we used the incompressibility constraint in order to 
compute F~ ; it states that det F = 1, so that 

r' = — 
Ar 



(2.5) 



5 



In our first deformation, the cylindrical tube is pressed into a cylindrical 
cavity with inner radius T\ > A and outer radius r 2 < B. It follows by 
integration of the equation above that 



r{R) = y ^ + a, (2.6) 
where now 

B 2 r 2 A2 2 B 2 -A 2 

B z — A 1 r 2 — r{ 

We compute the physical components of the left Cauchy-Green strain 
tensor B = FF l from (12.41) and find its first three principal invariants 1\ = 
tr B, I 2 = (det B)tr {B ), and I 3 = det B as 

I x = ( r ') 2 + {rg'f + (r/R) 2 + (rr) 2 + A 2 + {w'f, 

J 2 = (rA/i?) 2 + (ru/r - rt/A) 2 + (rw'/R) 2 + {R/rf + (1/A) 2 + (Rr/X) 2 , 

(2.8) 

and of course, J 3 = 1. 

For a general incompressible hyperelastic solid, the Cauchy stress tensor 
T is related to the strain through 

T = -pi + 2W X B -2W 2 B~ l , (2.9) 

where p is the Lagrange multiplier introduced by the incompressibility con- 
straint, W = W(Ii, h) is the strain energy density, and Wi = dW/dT. Hav- 
ing computed B 1 = (F t )~ 1 F~ 1 from (12.41) . we find that the components of 
T are 

T rr = -p + 2W 1 (r') 2 - 2W 2 [(r\/R) 2 + (rw'r - rg'X) 2 + (rw'/R) 2 ] , 

Tee = -V + 2IVi [(rg'f + (r/i?) 2 + (rr) 2 ] - 2H/ 2 (i?/r) 2 , 

T 2 , = -p + 2W 1 [\ 2 + (w'f] - 2W 2 [(1/A) 2 + (Rr/X) 2 ] , 

T re = 2W 1 (rr'g') - 2W 2 (w'r - g'\)R, 

T rz = 2W 1 (r'w > ) - 2W 2 [rRg'r - rRw'r 2 /X - rw'/(XR)] , 

T 6z = 2W 1 (rw'g' + rAr) + 2W 2 (r'Rr). (2.10) 

Finally the equilibrium equations, in the absence of body forces, are: 
div T = 0; for fields depending only on the radial coordinate as here, they 
reduce to 

dT rr T rr — Tae dT r e 2 m dT rz 1_ . . 

r L + ^r ee =Q _^L + _ Tr9 = Q ^ _^± + _ Trz = . 2.11 

dr r dr r dr r 
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2.2 Boundary conditions 

Now consider the inner face of the tube: we assume that it is subject to a 
vertical pull, 

T rz (A)=T A , T r9 (A) = 0, (2.12) 

say. Then we can integrate the second and third equations of equilibrium 
(12.1 |3 ; we find that 

T rz {r) = T j-T A 1 T r9 (r)=0. (2.13) 

The outer face of the tube (in contact with the glass in the cork/bottle 
problem) remains fixed, so that 

w(B) = Q, g(B) = Q, T rr (B) = T , (2.14) 

say. In addition to the axial traction applied on its inner face, the tube is 
subject to a resultant axial force N (say) and a resultant moment M (say), 

N= / T zz rdrd6, M= / T dz r 2 drd6. (2.15) 

JO Jri JO Jr\ 

Note that the traction T of ( 12. 14jl is not arbitrary but is instead deter- 
mined by the shrink fit pre- deformation ( 12.31) . by requiring that N = when 
T Q A = r = g = w = (This process is detailed in the next section for 
the neo-Hookean material.) Therefore T is connected with the stress field 
experienced by the cork when it is introduced in the bottleneck. 

In the rest of the paper we aim at presenting results in dimensionless 
form. To this end we normalize the strain-energy density W and the Cauchy 
stress tensor T with respect to /x, the infinitesimal shear modulus; hence we 
introduce W and T defined by 

W — T 

W = —, T=— . (2.16) 
Similarly we introduce the following non-dimensional variables, 



A^R_ri_w_a_ 
V = — ) H = — , Ti = — , w = — , a = — -, r = Bt, (2.17) 
' B B B B B 2 



so that r) < R < 1. Also, we find from (12. TH that 



1 — tj r% — r 
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Turning to our cork or shock absorber problems, we imagine that the inner 
metal cylinder is introduced into a pre-existing cylindrical cavity (this pre- 
caution ensures a one-to-one correspondence of the material points between 
the reference and the current configurations). In our upcoming numerical 
simulations, we take A = B/10 so that r] — 0.1; we consider that the outer 
radius is shrunk by 10%: r 2 = 0.9B, and that the inner radius has doubled: 
r\ = 2A; finally we apply a traction whose magnitude is half the infinitesimal 
shear modulus: |T A | = /i/2. This gives 

a~3.22xl0" 2 , A ~ 1.286, Tq — —0.5. (2.19) 

At this point it is possible to state clearly our main observation. A first 
glance at the boundary conditions, in particular at the requirements that 
g be zero on the outer face of the tube, gives the expectation that g = 
everywhere is a solution to our boundary value problem, at least for some 
simple forms of the constitutive equations. In what follows we find that for 
the so-called 'neo-Hookean' solids, g = is indeed a solution, whether there 
is a torsion r or not. However if the solid is not neo-Hookean, then it is 
necessary that g ^ when r ^ 0, and the picture becomes more complex. 
For this reason, we classify as 'purely academic' the question: which is the 
most general strain- energy density for which it is possible to solve the above 
boundary value problem with g = ? Indeed there is no 'real world' material 
whose behaviour is ever going to be described exactly by that strain-energy 
density (supposing it exists). Instead a more pertinent issue to raise for 'real 
word applications' is whether we are able to evaluate the importance of latent 
(secondary) stress fields, because they are bound to be woken up (triggered) 
by the deformation. 

3 neo-Hookean materials 

First we consider the special strain energy density which generates the class 
of neo-Hookean materials, namely 

W = (h- 3)/2, so that 2W 1 = 1, W 2 = 0. (3.1) 

Note that here and hereafter, we use the non-dimensional quantities intro- 
duced previously, from which we drop the overbar for convenience. Hence the 
components of the (non-dimensional) stress field in a neo-Hookean material 
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Figure 2: Pulling on the inside face of a neo-Hookean tube. Here the vertical 
axis is the symmetry axis of the tube. 



reduce to 

T rr = -p+ (r') 2 , Tee = ~P + (rgf + (r/R) 2 + (rr) 2 , 
T zz = -p + \ 2 + (V) 2 , T r0 = rr'g', 

T rz = r'w', T 6z = rg'w' + r\r. (3.2) 

Substituting into (I2.13P we find that 

w' = Xr \Tq / R } g' = 0, (3.3) 

and by integration, using f)2.14p . that 

w = ArxT A In R, g = 0. (3.4) 

In Figure 2a, we present a rectangle in the tube at rest. It is delimited 
by 0.1 < R < 1.0 and 0.0 < Z < 1.0. Then it is subject to the deformation 
corresponding to the numerical values of (12. 19ft . To generate Figure 2b, we 
computed the resulting shape for a neo-Hookean tube, using (12.21) . i \2.6L and 

flH3D. 

Now that we know the ful deformation field, see (12.21) and (13.41) . we can 
compute T rr — T 9 e from (13.21) and deduce T rr by integration of (I2.11l) i . with 
initial condition (12. 14ft ^. Then the other field quantities follow from the rest 
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of (13.21) . In the end we find in turn that 



T rr = — < In 



2A 



T 



Xr 2 R 2 
R 2 + a\ 
R 2 



+ {R 2 - I) 



a 



T rr + [ — + a) [ — + t 



A 



T rr + X [ 1 + 



R 2 
R 2 



r 2 (R 2 + aX) 
R 2 

X(R 2 + aX) 
R 2 



+ T 



(where we used the identity 1 + aX 



X(R 2 + aX) ' 
Arf, see (E6]) with R 



(3.5) 
1), and that 



T, 



i-6 



T 



n 



R 2 

-r- + a 
X 



--T, 



o J 



0Z 



>R 2 

Xt\/ — + a. 



(3.6) 



The constant T is fixed by the shrink fit pre- deformation (12.31) . imposing 
that N = when r = g = w = = 0, or 



(T + A 2 )(W 



Xr 2 R 2 



a(R 2 



2R 2 



l *^^ + r^ + :X) -Wf^x\ RdR = °- 

(3.7) 

Using this, and (12.151) . (13.51) . (I3.6p . we find the following expressions for the 
resultant moment, 

M = n{r\ - rf)Xr/2, (3.8) 



and for the axial force, 

N = 27rAr?|ln77| (T^f 



7T 

1' 



2\2 2 



(3.9) 



We now have a clear picture of the response of a neo-Hookean solid to the 
deformation (I2.2p . with the boundary conditions of §2 6. First we saw that 
here the contribution g(R) is not required for the azimuthal displacement, 
whether there is a torsion r or not. Also, if a moment M ^ is applied, then 
the tube suffers an amount of torsion r / proportional to M. On the other 
hand, if the tube is pulled by the application of an axial force only (N ^ 0) 
and no moment (M = 0), then r = and no azimuthal shear occurs at all. 

When we try to apply our results to the extrusion of a cork from the neck 
of a bottle, the following remarks seem to be relevant. From the elementary 
theory of Coulomb friction, it is known that the pulled cork starts to move 
when, in modulus, the friction force exerted on the neck surface is equal to 
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the normal force times the coefficient of static friction. In our case this means 
that 

\]\T rz {l)\ 2 + |T r ,(l)| 2 = f s \T„(1)\ = f s \T Q \ , (3.10) 

where fs is the coefficient of static friction. Using (12.121) and ( 12. 13ft . we find 
that the elements of the left handside of this inequality are 

T rz {l) = (n/r 2 )T^, T re {l) = 0. (3.11) 

Now, our main concern is to understand if it is better to apply a moment M ^ 
when uncorking a bottle, than to pull only. To address this question we 
note that the left handside of inequality (13.101) increases when \Tq-\ increases; 
on the other hand, combining (13. 8p and (13.91) . we have 

/(27rAr?|]n77|), (3.12) 

It is now clear, that for a fixed value of T^ 4 , in the case M ^ 0, it is necessary 
to apply an axial force whose intensity is less than the one in the case M — 0. 
Moreover, the equation above shows that (T^ 4 ) grows linearly with N but 
quadratically with M. With respect to efficient cork-pulling, the conclusion 
is that adding a twisting moment to a given pure axial force is more advan- 
tageous than solely increasing the vertical pull. Moreover, we observe that a 
moment is applied by using a lever and this is always more convenient from 
an energetic point of view. 

Recall that we made several simplifying assumptions to reach these re- 
sults: not only infinite axial length, incompressibility, and isotropy, but also 
the choice of a truly special strain energy density. In the next two sections 
we depart from the neo-Hookean model. 

4 Generalized neo-Hookean materials 

As a first broadening of the neo-Hookean strain-energy density (13. ip . we con- 
sider generalized neo-Hookean materials, for which the strain-energy density 
is a nonlinear function of the first invariant I\ only, 

W = W(h), (4.1) 

say. To gain access to the Cauchy stress components in this context, it 
suffices to take W 2 = and W\ = W in equations (12.101) . In particular, 
T r = 2rr'g'W, and the integrated equation of equilibrium (12.131) 9 gives 
g' = 0. Integrating, with (12.121) as an initial value, we find that 

9 = 0. (4.2) 



A\2 



N + 



vrA 2 



+ r 2 2 ) 2 



M 
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Hence, just as in the neo-Hookean case, azimuthal shear can be avoided 
altogether, whether there is a torsion r or not. We are left with an equation 
for the axial shear, namely ( 12. 131) i . which here reads 

2W'(h)w'(R) = ^-T A (4.3) 

Obviously the same steps as those taken for neo-Hookean solids may be 
followed here for any given strain energy density (14. ip . but now by resorting 
to a numerical treatment. Horgan & Saccomandi (2003a) show, through some 
specific examples of hardening generalized neo-Hookean solids, how rapidly 
involved the analysis becomes, even when there is only helical shear and no 
shrink fit. 

Instead we simply point out some striking differences between our present 
situation and the neo-Hookean case. We remark that I\ is of the form (I2.8l) i 
at g = that is, 

It follows that (14. 3p is a nonlinear differential equation for w', in contrast 
with the neo-Hookean case. Another contrast is that the axial shear w is now 
intimately coupled to the torsion parameter r, and that this dependence is 
a second-order effect (r appears above as t 2 ). 

A similar problem where the azimuthal shear has not been ignored, but 
the axial shear has been considered null i.e. w = has been recently consid- 
ered by Wineman (2005). 



5 Mooney— Rivlin materials 

In this section we specialize the general equations of §2 to the Mooney-Rivlin 
form of the strain-energy density, which in its non-dimensional form reads 

w= h-3 + nth -3) sothat 1 m 

2(1 + m) 1 + m 1 + m 

where m > is a material parameter, distinguishing the Mooney-Rivlin ma- 
terial from the neo-Hookean material (13.11) . and also allowing a dependence 
on the second principal strain invariant I2, in contrast to the generalized 
neo-Hookean solids of the previous section. 

Then the integrated equations of equilibrium (I2.13P read 

(R + mr 2 r 2 R + mr 2 / R) w' - (mT\r 2 R)g' = (1 + m)\r 1 T^, 
(mrA)w' - (1 + mAV = 0. (5.2) 
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First we ask ourselves if it is possible to avoid torsion during the pulling 
of the inner face. Taking r = above gives 



(R + mr 2 /R)w' = (1 + m)Ar 1 T ( f, g' = 0. (5.3) 

It follows that here it is indeed possible to solve our boundary value problem. 
We find 



A A(l + m) 1 

w = XriT n — - r In 

2(A + m) 



maX + (A + m)R 2 
maX + (A + m) 



g = 0. (5.4) 



However if r ^ 0, then it is necessary that g ^ 0, otherwise (15.21) 9 gives 
w' = while (I5.2l) i gives w' ^ 0, a contradiction. This constitutes the first 
departure from the neo-Hookean and generalized neo-Hookean behaviours: 
torsion (r ^ 0) is necessarily accompanied by azimuthal shear (g j^O). 

In the case r ^ 0, we introduce the function A = A(R) defined as 

A(R) = {R + mr 2 /R){\ + mX 2 ) + mr 2 r 2 R, (5.5) 

(recall that r = r(R) is given explicitly in (12. 6p .) We then solve the system 
(15.21) for w' and g' as 

w'=(l+ m)(l + m A 2 )A^rn, g' = m(l + m )A 2 -^|yrr 1 , (5.6) 

making clear the link between g and r. Thus for the Mooney-Rivlin mate- 
rial, the azimuthal shear g is a latent mode of deformation; it is woken up by 
any amount of torsion r. Recall that at first sight, the azimuthal shear com- 
ponent of the deformation (12. 2p seemed inessential to satisfy the boundary 
conditions, especially in view of the boundary condition g(l) = 0. However, 
a non-zero W2 term in the constitutive equation clearly couples the effects of 
a torsion and of an azimuthal shear, as displayed explicitly by the presence 
of r in the expression for g' above. 

It is perfectly possible to integrate equations (I5.6P in the general case, but 
to save space we do not reproduce the resulting long expressions. With them, 
we generated the deformation field picture of Figure 3a and Figure 3b. There 
we took the numerical values of ( 12.191) for a, A, Tq 4 ; we took a Mooney-Rivlin 
solid with m = 5.0; we imposed a torsion of amount r = 0.5; and we looked 
at the deformation field in the plane Z = 1 (reference configuration) and 
z = X (current configuration). 

Although the secondary fields appear to be slight in the picture, they are 
nonetheless truly present and cannot be neglected. To show this we consider 
a perturbation method to obtain simpler solutions and to understand the 
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Figure 3: Pulling on the inside face of a Mooney-Rivlin tube, with a clockwise 
torsion. 



effect of the coupling, by taking m small. Then integrating (15.61) . we find at 
first order that 

ZD I 

—j ~ (l+m)XlnR--m [t 2 R 2 + 2(1 + r 2 aX) In R - a\/R 2 - r 2 + aX] , 
riT 2 

5 



X 2 Tm\nR. (5.7) 



Hence, the secondary field g exists even for a nearly neo-Hookean solid (if m 
is small, then g of order m.) Interestingly we also note that the azimuthal 
shear g in (15. 7\\ varies in a homogeneous and linear manner with respect to 
the torsion parameter r and in a quadratic manner with respect to the axial 
stretch A, showing that that the presence of this secondary deformation field 
cannot be neglected when the effects of the prestress and of the torsion are 
both taken into account. To complete the picture, we use the first-order 
approximations 

2Wi ~ 1 - m, 2W 2 ~ m, (5.8) 
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, 2 



I 2 (r 2 r 2 r 2 X 2 I 

m 



( 1+2mA2 »7p-flbr + /p + fl-2i 



to obtain the stress field as 

T rr ~-p+(l-m) (r 1 ) 2 - m {(rX/R) 2 + [(rr) 2 + (r/R) 2 ] (XrjTf)* /R 2 } 
T ee ~ -p + (1 - m) [(r/R) 2 + (rr) 2 ] - m(R/r) 2 , 

T zz ~-p+(XT A n) 

+ (i _ m ) X 2 - m [(l/X) 2 + (Rt/X) 2 ] , 
T r e — rr' g' — mXrxT^r, 

T rz ~ (1 - m) (rV) + mAriT A [rRr 2 /X + r/(XR)] /R, 

T 9z ~ (1 - m) rAr + Xrr^g'/R + m(r'Rr). (5.9) 

Using this stress field it is straightforward, but long and cumbersome, 
to derive the analogue for a Mooney-Rivlin solid with a small m of relation 
(I3.12p (which was established for neo-Hookean solids.) However nothing 
truly new is gained from these complex formulas with respect to the simple 
neo-Hookean case, and we do not pursue this aspect any further. 



6 Conclusion 

In non-Newtonian fluid mechanics and in turbulence theory, the existence of 
shear-induced normal stresses on planes transverse to the direction of shear 
is at the root of some important phenomena occurring in the flow of fluid 
down pipes of non-circular cross section (Fosdick & Serrin 1973). In other 
words, pure parallel flows in tubes without axial symmetry are possible only 
when we consider the classical theory of Navier-Stokes equations or the linear 
theory of turbulence or tubes of circular cross section. 

In nonlinear elasticity theory, similar phenomena are reported. Hence 
Fosdick & Kao (1978) and Mollica & Rajagopal (1997) show that for gen- 
eral isotropic incompressible materials, an anti-plane shear deformation of a 
cylinder with non axial-symmetric cross section causes a secondary in-plane 
deformation field, because of normal stress differences. Horgan & Saccomandi 
(20036) give a detailed discussion of how the anti-plane shear deformation 
field couples with the in-plane deformation field in a generalized neo-Hookean 
solid. 

The appearance of what we called here latent deformations is quite gen- 
eral and common. For example it is known in compressible nonlinear elas- 
ticity that pure torsion is possible only in a special class of materials, but 
we know that torsion plus a radial displacement is possible in all compress- 
ible isotropic elastic materials (Polignone & Horgan 1991) (Here we signal 
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that 'possible in all materials' is not equivalent to 'universal', because the 
corresponding radial deformation differs from one material to another.) 

In this paper we give an example where axial symmetry holds, where the 
boundary conditions suggest that an axial shear deformation field is sufficient 
to solve the boundary value problem, and where nevertheless, the normal 
stress difference wakes up a latent azimuthal shear deformation. Moreover, 
because we are able to find some explicit exact solutions by some perturbation 
techniques, we are able to evaluate the importance of the latent deformation. 
Indeed, we show that if a certain constitutive parameter m (distinguishing 
a neo-Hookean solid from a Mooney-Rivlin solid) is zero or if the torsion 
parameter r is zero, then the solution to the boundary value problem can 
be found only in terms of the axial shear deformation field; if these two 
parameters are not zero — even if they are small — then the latent mode of 
deformation is quantitatively appreciable. 

In conclusion we suggest that it is not crucial to determine the class of 
materials for which a given deformation field is possible. Rather it is crucial to 
classify all the latent deformations associated with a given deformation field 
in such a way that this field is controllable for the entire class of materials. 
Indeed, no real material, even when we accept that its mechanical behaviour 
is purely elastic, is ever going to be described exactly by a special choice 
of strain-energy. Looking for special classes of materials for which special 
deformations fields are admissible may mislead us in our understanding of 
the nonlinear mechanical behaviour of materials. 

To finish the paper on a light note, we evoke a classic wine party dilemma: 
which kind of corkscrew system requires the least effort to uncork a bottle? 
Figure 4 sketches the two working principles commonly found in commercial 
corkscrews. The most common (on the left) relies on pulling only (directly 
or through levers) and the other type (on the right) relies on a combination 
of pulling and twisting. Notwithstanding the shortcomings of this paper's 
modelling with respect to an actual uncorking, the authors are confident 
that they have provided a scientific argument to those wine amateurs who 
favour the second type of corkscrews over the first. 
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